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FOREWORD 


The  development  of  symbolic  algebraic  manipulation  systems  for  use  on  the  digital 
computer^is  a  fairly  recent  technical  innovation  that  has  not  even  begun  to  be  fully  ex¬ 
ploited.  The  potential  for  application  of  these  systems  to  the  nonnumeric  resolution  of 
complex  problems  seems  boundless.  The  purpose  of  this  report  is  to  provide  the  results 
obtained  from  the  utilization  of  one  such  “state  of  the  art”^system,  MACSYMA,  to  provide 
high-order  expansion  formulae  for  the  classical  orbital  eccentricity  function  found  in  celestial 
mechanics  theory.  These  formulae  are  generally  applicable  to  low  eccentricity  orbits,  and 
the  relative  ease  with  which  they  were  generated  demonstrates  MACSYMA’s  analytical 
prowess.  This  report  was  reviewed  by  Dr.  R.  J.  Anderle  and  Mr.  R.  W.  Hill. 


Released  by: 


THOMAS  A.  CLARE,  Head 
Strategic  Systems  Department 


.  1  ,  ;  '  r  . 


flv jliabili >y  c 


y  cod »s 


NSWC  TR  84-355 


CONTENTS 

Page 

INTRODUCTION .  1 

THE  Ku  (e)  ECCENTRICITY  FUNCTION .  1 

8  pq 

DISCUSSION .  2 

REFERENCES .  3 

APPENDIX 

TABLE  OF  Kc  q(e)  ECCENTRICITY  FUNCTION  FORMULAE  FOR 

2  <  e  <  10,q  0  <  P  <  5,  AND  I  q  I  <  2 .  A-l 

DISTRIBUTION .  (1) 


NSWC  TR  84-355 


INTRODUCTION 


Although  digital  computers  have  been  performing  complex  numerical  computations 
for  several  decades,  their  application  to  the  solution  of  nonnumeric  problems  is  a  relatively 
recent  development.1,2  Several  sophisticated  symbolic  computer  algebra  systems  currently 
exist  and  can  serve  as  a  most  valuable  asset  to  the  scientific  researcher.3'4  Such  systems 
can  provide  general  solutions  to  a  wide  variety  of  problems  in  a  manner  that  eliminates  the 
need  for  tedious,  manually  performed  algebraic  operations,  while  simultaneously  minimiz¬ 
ing  the  likelihood  of  introducing  errors  into  the  solution. 

One  such  “state  of  the  art”  system  is  the  MACSYMA  computer  algebra  system  de¬ 
veloped  by  the  Mathlab  group  at  the  Massachusetts  Institute  of  Technology.5  In  this  work, 
MACSYMA  has  been  used  to  generate  expansion  formulae  for  the  eccentricity  function 
through  tenth  order  in  the  orbital  eccentricity.  These  formulae  are  presented  in  tabular 
form  for  2  <  &  <  10,  0  <  p  <  5,  and  |  q  j  <  2  and  may  be  used  to  verify  numerically 
generated  values  for  Kepq(e)  or  applied  directly  to  the  analytical  resolution  of  orbital  dy¬ 
namics  problems.  For  the  sake  of  completeness,  the  following  sections  also  include  a  mathe¬ 
matical  description  of  the  function  KCpq(e)  as  well  as  a  comparison  of  results  with  those 
obtained  by  other  workers. 


THE  K0  (e)  ECCENTRICITY  FUNCTION 

VI  p  q 


The  K  (e)  eccentricity  function  is  related  to  the  G^fe)  function  used  by  Kaula6 
through  the  expression 


G.  (e)  = 

t  ph 


'K*ph(c> 


(1) 


where.7  for  q  >  0: 


K,pq(e)=(-1)1  2k(  I+7) 
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DISCUSSION 


I  quations  2  through  4  of  the  previous  section  were  coded  in  the  MACSYMA  pro¬ 
gramming  language  and  processed  by  the  MACSYMA  system  to  generate  the  tables  of 
expansion  formulae  provided  in  the  attached  appendix.  Since  the  majority  of  satellite 
orbits  to  which  these  results  might  be  applied  are  of  low  eccentricity .  it  was  found  to  In¬ 
adequate  to  set  the  upper  limit  for  the  summation  index  k  in  liquations  2  anil  3  to 
5  and  constrain  the  q  subscript  to  i  q  |  <  2.M 
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The  MACSYMA  generated  results  were  compared  with  those  obtained  by  other  workers 
and  agreed  identically  with  those  fourth-order  formulae  presented  by  Giacaglia7  for  2  <  £ 
<  4,  0  <  p  <  4,  and  |  q  |  <  2.  The  MACSYMA  formulae  also  compared  identically  with 
Kaula’s6  fourth-order  expressions  for  G  (e)  for  2  <  C  <  4,  0  <  p  <  4,  and  |  q  i  <  2, 
except  for  the  cases  (C,  p.  q)  =  (4.  1.  -  1)  and  (4,  3.  1)  where  Kaula’s  second  term 

should  be  33c3  / 1 G  instead  of  33c2/  lb.  A  comparison  was  also  made  with  fourth-order 
results  for  G  (e)  for  0  <  p  <  2  and  jq|  <  2  given  by  Goldreich  and  Peale.9  The 
agreement  was  identical  except  for  the  cases  (C,  p,  q)  =  (2.  0,  -  1)  and  (2,  2.  1).  where 
their  first  term  should  be  -e/2  instead  of  e/2,  and  (2,  2,  -  1)  where  their  second  term 
should  be  -  1 23e3 / 1 6  instead  of  —  1 23e3 / 1  5- 

C'onsistency  checks  were  also  performed  upon  the  MACSYMA  formulae  using  the 
following  relations.10 

Kkpf+,  =  |  (38  -  4p  +  1)  +  0(e2 )  (7) 

and 

I 

K  =  -  (4p  -  8  +  1)  +  Ofe2 )  (8) 

General  overall  agreement  was  attained  to  this  order  of  accuracy. 
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APPENDIX  A 


TABLE  OF  KCpq(e)  ECCENTRICITY  FUNCTION* 
FORMULAE  FOR  2  <  £  <  10,  0  <  p  <  5, 
AND  |q|  <  2 
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